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We consider the existence of bulk chiral fermions around points of symmetry in the Brillouin zone 
of nonmagnetic 3D crystals with negligible spin-orbit interactions. We use group theory to show 
that this is possible, but only for a reduced number of space groups and points of symmetry that 
we tabulate. Moreover, we show that for a handful of space groups the existence of bulk chiral 
fermions is not only possible but unavoidable, irrespective of the concrete crystal structure. Thus 
our tables can be used to look for bulk chiral fermions in a specific class of systems, namely that 
of nonmagnetic 3D crystals with sufficiently weak spin-orbit coupling. We also discuss the effects 
of spin-orbit interactions and possible extensions of our approach to Weyl semimetals, crystals with 
magnetic order, and systems with Dirac points with pseudospin 1 and 3/2. A simple tight-binding 
model is used to illustrate some of the issues. 
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I. INTRODUCTION 

Electrons in the vicinity of the K points in graphene^] 
have linear dispersion relations and behave like mass- 
less chiral fermions. More concretely, the dynamics of 
electrons around these points is governed by the rel- 
ativistic, two-dimensional massless Dirac hamiltonian 



H tt 



CTyky, 



and many of the exotic electronic 



properties of graphene stem from this factP This also 
turns graphene into a potential laboratory for two- 
dimensional rclativistic dynamics, incorporating mass- 
less fermions, gauge fields and curved gravitational back- 
grounds.^ Moreover, optical lattices that could be used 
to simulate relativistic systems with trapped cold atoms 
can be fashioned after grap hene with control over the 
properties of the systemP-H Obviously three-dimensional 
analogs of graphene are potentially very interesting. 

Strictly speaking, the massless Dirac hamiltonian Hq 
describes only the low energy, orbital dynamics of elec- 
trons in graphene. As reviewed in Section III, spin-orbit 
coupling in graphene gives fermions a very small massP 
This mass is so small that for most practical purposes 
spin and orbital degrees of freedom decouple and Hq pro- 
vides an effective description of an enormous variety of 
phenomenaP 

In this paper we will consider 3D analogs of graphene, 
i.e., crystals with orbital electron dynamics governed 
by the 3D two-component massles Dirac hamiltonian 
Ho ~ va ■ k, also known as the Weyl hamiltonian. This 
hamiltonian describes massless chiral fermions, right- 
handed for v > 0, left-handed otherwise. Henceforth we 
will use the term 'orbital Weyl point' to refer to points 
around which the low energy dynamics in the absence of 
spin-orbit couplings is described by the 3D Weyl hamil- 
tonian, with an additional twofold degeneray due to elec- 
tron spin. We will also speak of 'bulk chiral fermions', 
keeping in mind that, as in grapheneP they are exactly 
chiral in the limit of vanishing spin-orbit interactions. 
Also note that it is pseudospin, not electron spin, which 



is parallel (or antiparallel) to k in the chiral limit and 
that, in these systems, pseudospin is purely orbital in 
origin. 

These should be distinguished from other systems 
where the total hamiltonian, including electron spin and 
spin-orbit interactions, adopts the form of the Weyl 
hamiltonian. These include surface states in topologi- 
cal insulators^ as well as novel three-dimensional 'Weyl 
semimetals'J23E^ The spectrum of these systems, unlike 
graphene and its 3D analogs considered in this paper, re- 
mains gapless for arbitrary values of the spin-orbit cou- 
plings. Actually, some Weyl semimetals have strong spin- 
orbit interactions.^ Possible extensions of our methods 
to Weyl semimetals will be considered in the last Section. 

It is well known that Weyl points have topological 
propertiesp^Hm and no fine-tunning or symmetries are 
usually necessary for their existence. But symmetry, 
while not necessary, can sometimes be sufficient for the 
existence of Weyl points. That is what we show in this 
paper, where we investigate the role played by the space 
groups of crystals with time reversal symmetry (TRS) in 
the existence of orbital Weyl points. 

The main results of this paper are summarized in Ta- 
bles HJn] Only crystals with one of the 19 space groups 
in these tables can have orbital Weyl points at points of 
symmetry. Moreover, crystals with space groups in Ta- 
ble |l] must have orbital Weyl points at the listed points, 
irrespective of the actual crystal structure. For crystals 
with space groups in Tabic [j*T| the situation is only slightly 
different: At the listed points both orbital Weyl points 
and non-degenerate bands with quadratic dispersion re- 
lations are possible. These results are relevant to non- 
magnetic 3D crystals with sufficiently weak spin-orbit 
interactions and to cold atoms in optical lattices. 

The rest of the paper is organized as follows. Our main 
results, contained in Eq. ^ and Tables IJlI are explained 
in Section II. Section III considers the effects of spin- 
orbit interactions on the orbital Weyl points, and a sim- 
ple tight-binding model is constructed and analyzed in 



Section IV. Possible extensions of our approach to Weyl 
semimetals, crystals with magnetic order and other types 
of Dirac points are considered in Section V. An outline 
of the methods used to obtain Tables |I|TT| is given in the 
Appendix. 

II. ORBITAL WEYL POINTS AND CRYSTAL 
SYMMETRY 

Our strategy is based on the fact that the form of the 
hamiltonian in the vicinity of a point of symmetry K\ is 
strongly constrained by the symmetries of the poi nt in 
questionPHU These include Q Kl —the little grourJMSil 

of the vector K\ — and combinations of TRS with space 
group elements. Since we are interested in systems with 
very weak spin-orbit interactions, we will consider first 
the structure of the orbital or spin-independent part of 
the hamiltonian. The transformation properties of or- 
bital wavefunctions are described by single- valuecP^ rep- 
resentations of the space group. 

As explained in the Appendix, we have carried out a 
survey of all the single- valued irreducible representations 
of the 230 space groups at points of symmetry in the Bril- 
louin zone (BZ). We find that, for most space groups, 
the constraints on the form of the hamiltonian around 
points of symmetry are incompatible with the structure 
of the Weyl hamiltonian. The comparatively few ex- 
ceptions are listed in Tables [T||TT] In all cases, two elec- 
tronic bands transforming according to a single-valued 
irreducible representation (IR) of Qk ± are degenerate at 
the point of symmetry K\. Near the point of symme- 
try, i.e., for K = Ki + k, the degeneracy is broken by 
fc-dependent terms and the hamiltonian takes the form 

H(k) = v x cr x k x + v y a y k y + v z a z k z + 0(k 2 ) (1) 

where v x = v y for uniaxial crystals and v x = v y = v z 
for cubic crystals. After appropriate rescalings of the 
components of k for non-isotropic crystals, this is just 
the Weyl hamiltonian H ~ va ■ k, with the sign of v 
equal to the product of the signs of Uj. The points of 
symmetry and IRs where this happens are listed in the 
last column of Tables ITlfTTl in standard notationPSl 

TRS reverses the sign of K. In those cases where — K\ 
is not equivalent to K\ , we get a copy of the Weyl hamil- 
tonian at the mirror point —Ki and fermions have, in 
addition to the pseudospin index associated to the Pauli 
matrices o~i , a 'valley' index. As shown in the Appendix, 
the total hamiltonian is then given by the 4x4 matrix 

This describes two degenerate massless fermions of the 
same chirality, right-handed for v > 0, left-handed oth- 
erwise. Somewhat surprisingly, we find that this doubling 
continues to take place even when — K\ = K\, i.e., for 
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Space Group 




ITVS 


214 


74i32 


O 8 


Pl,P2,P 3 


213 


P4i32 


Cf 


Rl , R 2 , R 3 


212 


P4 3 32 


o 6 


Rl , i?2 , R 3 


199 


72i3 


T 'l 


Pl,p 2 ,p 3 


198 


P2i3 




Ri,R 2 ,R 3 


098 


J4i22 


,-,1(1 


Pi 


096 


P4 3 2!2 




A 1 ,A 2 


092 


P4i2i2 


Di 


A X ,A 2 


024 


72i2i2i 




Wi 


019 


P2i2i2i 


Di 


Ri 



TABLE I: Cubic, tetragonal and orthorhombic space groups 
with orbital Weyl points. The small IRs are all 2d (except for 
i? 3 of 212 and 213, which is Ad) and refer to the symmetry 
points P(i, |, |), R(l §, |), A{\, |, §) and W{\, \, \), with 
components in the conventional basis of RefP^The stars have 
two vectors (Ki,— Ki) for body-centered (7) lattices, and a 
single vector Ki = —K\ for simple (P) lattices. 

TRS invariant momenta. In that case, Eq. ([2| describes 
two distinct fermions of the same chirality at the same 
point in the BZ. There is still a 'valley' index but, unlike 
in graphene, it can not be associated with two different 
points in the BZ. This happens for the six space groups 
with simple (P) Bravais lattices in Table [i] 

The groups in Table [I] have one important feature in 
common: The IRs in the last column of the table in- 
clude all the IRs at the point of symmetry. This means 
that, at that point, all the bands must form degener- 
ate pairs with Weyl hamiltonians. In other words, any 
crystal with space group in Table |T] will have bulk chiral 
fermions described by Eq. ^ , irrespective of the concrete 
crystal structure. On the other hand, the space groups 
in Table [TT| have, besides the listed small IRs K s and iJ 3 , 
which are two-dimensional and give rise to orbital Weyl 
points, other one-dimensional small IRs (K\,K2, H\,H-i) 
not related to Weyl points. In this case, both orbital 
Weyl points and non-degenerate bands with quadratic 
dispersion relations are possible at the listed points of 
symmetry. 
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182 
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Dl 


K 3 


181 


P6 4 22 


Dl 


K 3 ,H 3 


180 


P6 2 22 


Dl 


K 3 ,H 3 


179 


P6 5 22 


Dl 


K 3 


178 


P6i22 


Dl 


K 3 


177 


P622 


Dl 


K 3 ,H 3 


154 


P3 2 21 


Dl 


K 3 ,H 3 


152 


P3i21 


Di 


K 3 ,H 3 


150 


P321 


Di 


K 3 ,H 3 



TABLE II: Hexagonal and trigonal space groups with orbital 
Weyl points. The small IRs listed are all 2d and refer to the 
symmetry points K(^, |, 0) and 77(|, § , |), with components 
in the conventional basis of Ref The stars have two vectors 
(Ki, —Ki) in all cases. 
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FIG. 1: Effects of spin-orbit coupling on the orbital Weyl 
point of a cubic crystal (top) and graphene (bottom) in arbi- 
trary units. 



We close this section by pointing out some special fea- 
tures in Tables jlflll The first one is that all th e groups 
m Table [T] are subgroups of the first entry J2MZI t ne cubic 
space group 214 /4i32 and, despite the use of differ- 
ent conventional names (P, R, A, W), they also share the 
point of symmetry. Indeed, the cartesian coordinates for 
all the points in Table [I] can be written as 



6' c' 



(3) 



in terms of the unit cell constants, with b = c for uniaxial 
crystals and a = b = c for cubic crystals. As a result, 
one can begin with any 3D lattice with space group 214 
and reproduce all the cases in Table [l] by suitable defor- 
mations. The second somewhat surprising feature is that 
all the stars have just one or two vectors. As a conse- 
quence, these crystals have only one or two orbital Weyl 
points degenerate in energy. This should be contrasted, 
for instance, with the case studied in Ref.,^ where 24 
Weyl points (away from points of symmetry) are present 
at the Fermi energy. 



III. SPIN-ORBIT INTERACTIONS 

Strictly speaking, our analysis so far applies only to 
'spinless electrons'. It is well known that spin-orbit inter- 



actions in graphene open a gap and give fermionic excita- 
tions a small massPln the case of graphene, the intrinsic 
spin-orbit hamiltonian is proportional to a z <g> s z ® t z , 
where s z and t z are Pauli matrices for electron spin 
and valley indices respectively. Around each valley, the 
hamiltonian can be written in terms of 4 x 4 matrices 



H so - 



-v{a x k x 



0A 



(4) 



where on = Ui ® l s , j3 = ±<r z <E) s z and A is the strength of 
the spin-orbit coupling. The matrices satisfy the Clifford 
algebra 



{a l ,<x,} = 2%, {a i; /3} = 0, 



(5) 



This identifies Eq. Q as the Dirac hamiltonian for 
4-component fermions with mass m — A and spectrum 
E ± =±VA 2 + v 2 k 2 . 



For the space groups in Tables [ipi| the most general 
fc-independent spin-orbit hamiltonian compatible with 
spatial symmetries and TRS takes the form 



H so = - (A x a x 



-AyOy 



s y + A z a z 



s z )®l T (6) 



where — A y for uniaxial crystals and A^ = A y = A z 
for cubic crystals. Actually, this valley independent 
form of the spin-orbit interaction is valid on the basis 
(ex, 62, ie*}, — ie*) of orbital wavefunctions introduced in 
the Appendix. If one uses the more conventional basis 
{e\, e2, ej, ej), then one has to append the valley matrix 
t z to the x and z components in Eq. j6j). Note that, in 
the conventional basis, instead of Eq. (pF, we would have 
Eq. ( [A6| . 

The spectrum of the total hamiltonian va ■ k + H so can 
be computed numerically and one finds that, in general, 
gaps are generated and all fermionic excitations acquire 
masses. Cubic crystals, where H so is isotropic and de- 
pends on a single parameter A, are an exception and can 
be treated analytically. In this case the spectrum is given 
by 



El 



± vk 



(7) 



and contains massless excitations. This spectrum is rad- 
ically different from that of the Dirac hamiltonian ap- 
propriate for 2D graphene. Note, in particular, that the 
linear bands E±, together with ET, form a triplet (see 
Fig. [T]), following the usual rules for addition of angu- 
lar momenta with L = S = 1/2 and J = L + S. This is 
only natural: assembling pseudospin and spin into a 4- 
componcnt object is equivalent to taking the Kronecker 
product of two j = 1/2 irreducible representations of the 
50(3) rotation group and this product decomposes ac- 
cording the rules of angular momentum compositionP^l 
Thus, unlike in 2D where strong spin-orbit interactions 
simply destroy the orbital Weyl points and turn massless 
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fermions into massive excitations, here we also get Dirac 
points with J = 1. In the next Section we will present 
a model that in the absence of spin-orbit has, besides 
orbital Weyl points, Dirac points with pseudospin-one. 
The points with pseudospin-one would be split by strong 
spin-orbit interactions into Weyl points with J = 1/2 
an d nov el Dirac points with J = 3/2. Subduction of 
IRs 24 25 can be used in principle to extend the analysis 
to non-cubic groups. 

Henceforth we will assume that spin-orbit interactions 
are weak and can be ignored. In this limit the dynamics is 
well described by the 3D Weyl hamiltonian -albeit with 
an additional two-fold degeneracy due to spin- and, as 
a consequence, the systems considered in this paper may 
share some of the properties of Weyl semimetals. 

IV. A TIGHT-BINDING EXAMPLE 

As a practical application, we present a tight-binding 
model with space group 214 74i32(0 8 ), the first entry in 



Table [TJ According to our previous analysis, any lattice 
with this space group must have orbital Weyl points at 
P. Here we consider a lattice with four atoms per primi- 
tive unit cell, with cartesian coordinates fx — a/8(l, 1, 1), 
r 2 = a/8(l,l,3), f 3 = a/8(3,l,5) and r 4 = a/8(3, 1, 7)0 

To each atom we associate a Bloch function 



= e* (?I+ '~V(r - n - i) (8) 



where the sum runs over all the points in the Bravais lat- 
tice and (p(r) is an s-wave atomic orbital. Each atom has 
three nearest neighbors (NN) , with bonds parallel to the 
three cartesian planes. In terms of the reduced cartesian 
components of the wave vector k — 2tt /a(k Xl k yi k z ) the 
NN tight-binding hamiltonian is given by 



H(k) = t 
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FIG. 2: Bands for the cubic lattice in the NN approximation 
for t = — 1. The BCC Brillouin zone with its points and lines 
of symmetry can be seen in Fig. [3] 

where t < is the hopping parameter and the diagonal 
on-site energy has been set to zero. The hamiltonian can 
be diagonalized numerically and the resulting bands are 
shown in Fig. [2j 

The existence of linear bands at point P is obvious in 
Fig. [2j We can use standard group theory techniques 



to confirm that they actually are Weyl points, as pre- 
dicted. The four Bloch functions 4>i form the basis of 
a reducible representation % that can be decomposed 
into small IRs as . . . , $4) = Pz(ei, e2)+Ps(ui, U2). 

Up to normalizations, the symmetry- adapted vectors are 
given by e x ~ (1, j3 2 , -1, i/? 2 ), e 2 ~ (-fa, -1, -#2, i), 




FIG. 3: Brillouin zone for BCC crystals.^ 21221 
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with /?2 = (1 + \/3)(l —i)/2, and identical expressions for 
Ul, «2 with 02 replaced by ^3 = (1 — v3)(l — i)/2. Using 
a unitary transformation Up to change to the symmetry 
adapted basis and expanding around the point P yields 



U P H(k)U F 



irt 



f+CT-fc H 23 (k) 



where H 2 z{k) is given by 



2 I w'fcj + iujky 



uj*k x — iuk 
-k. 



0(k 2 ) 
(10) 

(11) 



with id — e 27ri / 3 . There are thus two orbital Weyl points 
at P with different energies ±tv3 and opposite chirali- 
ties. The existence of orbital Weyl points of opposite chi- 
rality is of course to be expected from fermion doubling, 28 
which requires the net chirality of the BZ to vanish, al- 
though the fact that they appear at coincident points is 
peculiar to this model. Note also the linear couplings be- 
tween the two points. Due to the split in energies, these 
couplings contribute corrections 0(k 2 ) to the 2 x 2 ef- 
fective hamiltonians around the orbital Weyl points and 
do not spoil their structure. A similar expansion around 
—K\ confirms that, due to TRS, each orbital Weyl point 
is degenerate in energy with another point of the same 
chirality. Restoring the lattice constant yields a Fermi 
velocity vp = ^7=- Note however that, as no symmetry 
connects the orbital Weyl points with different chirali- 
ties and energies, going beyond the NN approximation is 
expected to give different Fermi velocities for them. 

Fig. [2] exhibits linear bands around the T and H points 
as well. Their nature is, however, very different from that 
of the orbital Weyl points at P. Let's consider, for the 
sake of concreteness, the T point. In this case, the rep- 
resentation associated with the Bloch functions decom- 
pose into IRs of dimension one and three, V. = Ai + T 2 . 
Transforming to the appropriate symmetry-adapted ba- 
sis and linearizing yields 



UlH{k)U T = t 



( 3 
-1 

iirk z 

\ —iirky 











iirk z inky 
—iirkr- 



-1 
ink- 



-1 



0(k 2 ) 



(12) 



Up to a constant energy, the 3x3 block can be written 



H T2 (k) =irtJ-k 



(13) 



where (Ji)jk — —i^ijk are spin-1 matrices. Thus, around 
this Dirac point, electrons behave more like massless 
spin-one particles, with spectrum E(k) — 0,±vp\k\ and 
Vp = \a\t\. Indeed, one can check that, while the E = 
component is longitudinally polarized, the other two are 
transverse, just like the propagating components of a 
photon. Pseu dosp in-one Dirac points have been reported 
in some twc^HH] anc j three-dimensional systems. 

A look at Fig. [2] shows that, even if the Fermi level 
coincides with one of the orbital Weyl points at P, band 
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FIG. 4: Bands for the cubic lattice with a tetragonal distor- 
tion (e = 0.4) for t — —1.25. See the main text for details. 



overlap will cause the dynamics to be dominated by large 
electron (or hole) pockets. The existence of band over- 
lap can be traced in this case to the 3-fold degeneracies 
in Fig. [2j which force one of the bands arising from the 
orbital Weyl points to bend over. As 3d IRs exist only 
for cubic groups, we may modify the model by making 
the hopping parameters ij_ associated with bonds paral- 
lel to the OXY-plane different from the rest, t± = et. 
This reduces the symmetry to the tetragonal subgroup 
98 I4i22 and eliminates the 3-fold degeneracies at T and 
H. Fig. [4] shows the bands for e = 0.4. In the absence 
of spin-orbit interactions the system would behave like a 
gapless semiconductor with massless carriers for 1/4 and 
3/4 fillings, with positive or negative chirality depending 
on the filling. 



V. DISCUSSION 

In this paper we have studied the interplay be- 
tween crystal symmetry and the existence of bulk chi- 
ral fermions in nonmagnetic 3D crystals with weak spin- 
orbit coupling. We have shown that the space group 
plays a determinant role, as summarized in Tables [T|[TT] As 
Weyl semimetals depend on the existence of magnetic or- 
der or strong spin-orbit interactions,^ we have explored 
an entirely different corner in the space of candidate 3D 
crystals with bulk chiral fermions. We have also consid- 
ered the effects of spin-orbit interactions and shown that, 
unlike in 2D, these may give rise to new critical points 
supporting massless fermions for different values of the 
pseudospin. For sufficiently weak spin-orbit interactions 
the dynamics is well described by the 3D Weyl hamilto- 
nian and the systems considered in this paper may share 
some of the properties of Weyl semimetals. 

There are two obvious uses for the information in Ta- 
bles p|ll| First, one can look for orbital Weyl points in 
nonmagnetic crystals with weak spin-orbit interactions 
and space groups in the tables, either in theoretically 
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computed electronic bands or experimentally. The other 
use is the design of 3D lattices with Weyl points. This 
may allow for a physical realization of massless chiral 
fermions with cold atoms in optical lattices. Note that 
the physical relevance as well as the feasibility of detect- 
ing chiral fermions in actual crystals will be affected by 
structure dependent features such as the position of the 
Fermi level and the amount of band overlap. But know- 
ing that the bulk chiral fermions have to be there is ob- 
viously a good starting point. As shown in the example 
of Section IV, we can then try to engineer the required 
properties by modifying the initial system. 

There are several possible extensions to this work. One 
is to consider the existence of orbital Weyl points away 
from points of symmetry. All the points in a line of sym- 
metry, except for the endpoints, have the same group 
Qk i • Therefore, symmetry alone can not imply the exis- 
tence of bulk chiral fermions at a particular point in the 
line, although it will indicate whether this is possible at 
all. But it can, in some cases, imply the existence along 
the line of 'semi-Dirac' points, where the dispersion re- 
lations are linear only for some directions. For generic 
points in the BZ, symmetry alone has l ittle t o say and a 
different kind of analysis may be usefulP^EH 

Sometimes not two, but three bands become degener- 
ate at a Dirac point PSl This is the case with the T and 
H points in the example of Section IV. This is possi- 
ble for some cubic space groups, and group theory can 
be used to determine the IRs and points of symmetry 
where this may happen. As shown in Section III, turn- 
ing on spin-orbit interactions will give rise to novel Dirac 
points with pseudospin 3/2. Other, more complicate lin- 
ear hamiltoniansSSHSZl a re also possible and may be phys- 
ically relevant. Group theory can be used to classify or 
even predict the existence of the different types of Dirac 
points. 

In this paper we have analyzed all the single- valued ir- 
reducible representations at points of symmetry in the 
BZ. These are appropriate for orbital degrees of free- 
dom. By considering instead double-valuecP^ IRs we 
could study the existence of Weyl points where the total 
hamiltonian, including electron spin and spin-orbit in- 
teractions, takes the form of the 3D Weyl hamiltonian. 
Thus, we could extend our analysis to TRS invariant 
Weyl semimetals. Double- valued IRs have been recently 
applied to the study of 'Dirac semimetals'P^ 

So far we have restricted ourselves to TRS invariant 
crystals. The reasons are mostly practical. The sym- 
metries of crystals with magnetic order are classified by 
the 1651 magnetic space groups^, instead of the 230 
ordinary (Fedorov) space groups that classify crystals 
with TRS. The amount of work required to examine the 
points of symmetry and irreducible corepresentations^ 
for all the magnetic groups is, obviously, much greater. 
Moreover, unlike ordinary space groups, there are few 
databases with the magnetic space groups of crystals with 
magnetic order. 

Nevertheless, some of the results in this paper can also 



be applied to spinless electrons in crystals with mag- 
netic order. The reason is that, by construction, the 19 
entries in Tables HlHIl describe situations where the 3D 
Weyl hamiltonian is invariant under the 'grey' or 'type 
IF Shubnikov space group^ associated to an ordinary 
(Fedorov) space group by the addition of the TRS oper- 
ation. Now, all the magnetic space groups deriyed from 
the Fedorov space group with the BNS settings^EH are 
subgroups of the grey group. As a consequence, the cor- 
responding Weyl hamiltonian is automatically invariant 
under all the magnetic groups derived from the ordinary 
space groups listed in our tables. For instance, the Weyl 
hamiltonians at the P point of group 214 are automati- 
cally invariant under the derived magnetic groups 214.68 
and 214.69 (in the BNS settings). 

One still has to check that the degeneracy of the two 
bands at the point of symmetry, necessary for the exis- 
tence of the Weyl point, is maintained under the lower 
symmetry of the magnetic space group. If this is not the 
case but the effects of the magnetic order are small, the 
Weyl point will survive, but move away from the point of 
symmetry. The main difference when dealing with crys- 
tals with magnetic order is that we can not exclude the 
possibility of finding bulk chiral fermions around points 
of symmetry for space groups not listed in Tables IJlI 



That happens whenever the Weyl hamiltonian is invari- 
ant under the magnetic subgroup of the grey space group, 
but not under the grey space group itself. 
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Appendix A 

In this Appendix we outline the methods used to ob- 
tain Tables [TpT] and Eq. ([2| . Symmetry operations be- 
longing to the space group Q will be written g = {a\v\, 
where a and v denote the rotation and translation parts 
respectively.^ Let {ei(K)} be a basis of orbital wave 
functions that transform linearly under the action of Q, 
ei(K) — ¥ ej(aK)Rji(g), where we sum over repeated in- 
dices. Invariance of the hamiltonian under Q means that, 



H)i/, , where ip g is the 



for any wavefunction ip, (H)tp 
transformed of ip by the group element g. Expanding this 
condition on the basis {ei(K )} yields, in matrix notation 



R\g)H{aK)R(g) = H{K) 



(Al) 
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where Hij(K) = (e{(K )\H\ej(K)). Time reversal 9 is an 
antiunitary operation that acts on orbital wavefunctions 
by complex conjugation, ej(K ) — » e .j{K)* — ej(—K)Qji, 
where Qij is a unitary matrixp2!2SI an d reverses the mo- 
mentum K . Invariance under 9 implies 



Q ] H{-K)Q = H*{K) 



(A2) 



We will also have to consider combined antiunitary op- 
erations of the form 9g. In this case invariance of the 
hamiltonian implies 



T^g)H(-aK)T(g) = H*(K) 



(A3) 



where T(g) = QR*(g). Eqs. ( |A1|A2|A3[ ) become pow- 
erful constraints on the form of the hamiltonian when 
we take K = K x + k in the neighborhood of a point of 
symmetry K\ and consider a power expansion in k. 



Consider for instance Eq. ( Al I with g restricted to the 
little group of K\, Gku i-e. 



RHg)H(K 1 



aKi 

ak)R(g) = H{K X 



k) 



By definition, the matrix a belongs to the vector repre- 
sentation y) 23 l 24 llf the basis functions {ei(K)} belong to 
the small IR R of Qk 1 , terms of order n in k in an expan- 



sion of the l.h.s. of Eq. (A4| will transform according to 
the product R* x Rx [V ™, where R* is the complex con- 
jugate of R and [V] n d enotes the n-th symmetric power 
of the representation y.HSM] Then one can use standard 
group theory techniques to determine, order by order in 
k, the most general form of the hamiltonian compatible 
with the symmetries of the vector K\ . 

In particular, a necessary condition for the existence 
of the Weyl hamiltonian, which is linear in k, is that 
the vector representation V is contained in the product 
R* x R for some 2d IBP R of g Kl . We have checked 
this condition on all the single-valued 2d and Ad IRs at 
the points of symmetry of the 230 space groups to obtain 
a first list of candidatesJS Polar groups have been dis- 
carded from the outset, as they couple one component of 
the momentum to the unit matrix and the result is in- 
compatible with the structure of the Weyl hamiltonian. 
The amount of work involved at this stage has been sub- 
stantially reduced thanks to the use of 'abstract groups' 
in RefpSl In essence, the abstract groups represent classes 
of isomorphic little groups Qk x and their number is much 
smaller than the number of little groups. 

In the next step each candidate IR has been checked for 
invariance of the corresponding hamiltonian under TRS. 



This inv olves using Eq. (A2| whenever —K\ = K\ and 
Eq. (A3) if K\ is not equivalent to — K\ but there exists 
a space group element g — \a\v\ € Q such that —aKi = 
K\. The IRs that pass this last test are listed in the 
last column of Tables [TUTT] These IRs have an important 
property: There is always a basis where the matrices R(g) 
coincide, up to ^-dependent phases, with the j = 1/2 
rotation matrices 



i?i/ 2 (n</>) = exp I - 



ia ■ h 4> 



(A5) 



where <fr is the angle around the unit vector n. As a con- 
sequence, we can always transform to a basis where, up 
to appropriate rescalings of the components of k for non- 
isotropic crystals (see Eq. ([!])), the linear hamiltonian 
takes the standard Weyl form H(K\ + k) ~ v a ■ k. 

When —K\ is not equivalent to K x , the basis at these 
two points can be related by TRS. Choosing as basis 
at —Ki the complex conjugate of the one at K%, i.e., 



(A4) taking as 4d basis (ei, e2, e*, ejj), and using Eq. (A2) 



gives H{—Ki + k) ~ —va* ■ k. The off-diagonal blocks 
between Ki and —K\ vanish by translation invariance, 
and we have 



H(k) 



a ■ k 








+ 0{k 2 ) (A6) 



We can make the symmetry between Ki and —K\ ob- 
vious by using the SU(2) transformation a y a*a y = —a 
to change the basis at —K\ so that the Ad basis is 
(ex, e2,ie2, — On this basis the hamiltonian takes 

the form given in Eq. 

The case — K x = Ki is more subtle, and we have two 
possibilities. If R is a real 2d IR, then TRS applies 
(ei,e2) onto itself, with e* = ej. In this case, = 1 



and Eq. (A2) requires H(—k) = H(k)* 1 which is not sat- 
isfied by the Weyl hamiltonian. Thus 2d real IRs are 
excluded from Table [I] For complex and pseudoreal 2d 
IRs and for the real Ad R 3 of 212 and 213, (e*,el) are 
linearly independent of (ei, e 2 ). We can take as Ad basis 
either (ei, e2, e*, e|) or (ei, e 2 , ie\, — ie*), and everything 
proceeds as in the previous case. The main difference 
is that now the off-diagonal blocks need not vanish by 
translation symmetry. However, a detailed case by case 
analysis using space group and TRS invariance shows 
that the off-diagonal terms are at least of 0(k 2 ). This 
completes the proof of Eq. ^ . 
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